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ABSTRACT. The third symposium on Remote Sensing of Snow and Ice 
organized by the International Glaciological Society, took place in Boulder' 

Colorado, 17-22 May 1992. As part of this meeting a total of 21 papers was 
presented on snow and ice applications of Advanced Very High Resolution 
Radiometer (AVHRR) satellite data in polar regions. Also during this meeting a 
NASA sponsored Workshop was held to review the status of polar surface 
measurements from AVHRR. In the following we have summarized the ideas and 
recommendations from the workshop, and the conclusions of relevant papers given 
during the regular symposium sessions. The seven topics discussed include cloud 
masking, ice surface temperature, narrow-band albedo, ice concentration lead 
statistics, sea-ice motion and ice-sheet studies with specifics on applications 
algorithms and accuracy, following recommendations for future improvements. In 
general, we can affirm tbe strong potential of AVHRR for studying sea ice and snow 
covered surfaces, and we highly recommend this satellite data set for long-term 
monitoring of polar process studies. However, progress is needed to reduce the 
uncertainty of the retrieved parameters for all of the above mentioned topics to make 
this data set useful for direct climate applications such as heat balance studies and 
others. Further, the acquisition and processing of polar AVHRR data must become 
better coordinated between receiving stations, data centers and funding agencies to 
guarantee a long-term commitment to the collection and distribution of high quality 
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1. INTRODUCTION 

The Advanced Very-High Resolution Radiometer 
(AVHRR) is a satellite sensor with strong potential for 
studying sea ice and snow covered ice surfaces. Its visible 
and infrared channels, listed in Table 1, offer data with 
direct applicability to problems of heat balance. Its 
spatial resolution of 1 km at nadir elucidates the complex 
structures of leads and floes of sea ice and flowbands of ice 
sheets. Its frequent polar coverage and broad swath are 
sufficient to resolve synoptic and even diurnal changes of 
surface temperature and lead patterns, and to track sea- 
ice motion. The decades-long availability and planned 
continuation of this sensor make it a primary source of 
data that can be use to monitor global climate. 

Table 1. A VHRR wavelengths in microns. Swath width is 
2580 km High Resolution Picture Transmission 
(HRPT). Ground resolution is 1.1km at nadir for 
Local Area Coverage (LAC) and HRPT. Global Area 
Coverage (GAC) data have a ground resolution of 4 km 
(after Kidwell, 1991) 


Band Spectral region 

Wavelength 

Wavelength 


NOAA-6,8, 10 

NOAA-7,9, 11, 12 



fim 

/im 

1 

Visible 

0.58-0.68 

0.58-0.68 

2 

Near infrared 

0.725-1.10 

0.725-1.10 

3 

Infrared 

3.55-3.93 

3.55-3.93 

4 

Thermal infrared 

10.50-11.50 

10.30-11.30 

5 

Thermal infrared 

10.50-11.50 

11.50-12.50 


Despite this promise, AVHRR has not been broadly 
applied to polar scientific problems. A number of studies 
have been undertaken to show how geophysically 
important surface parameters can be inferred using 
AVHRR data, beginning two decades ago with three- 
day minimum-brightness estimates of ice concentrations. 
Most studies, however, have worked in a demonstration 
and/or case study mode, and have not been developed 
into a monitoring capability. In contrast, consider the 
widely accepted data set to study the cloud cover 
(International Satellite Cloud Climatology Project, 
ISCCP), the variability of biomass through a vegetation 
index, and the sea surface temperature. There is no 
analogy to these data sets for any polar surface variable 
derived from AVHRR data. A reason for this lack of 
applications to the polar regions has been and continues 
to be cloudiness that obscures the surface for most of the 
summer and much of the winter. Although there are holes 
in the cloud cover through which the surface can be 
viewed, differentiating cloud from surface over snow and 
ice is far more difficult than at lower latitudes where the 
albedo and temperature of land and water are quite 
different from those of clouds. Another problem is the 
diverse number of institutions and agencies where 
AVHRR data are archived, in multiple formats and 


media, and for various periods, all without inter-related 
catalogs. For this reason, the utilization of polar AVHRR 
data is poor. 

There is some feeling that the polar community needs 
better access to polar-gridded, calibrated, cloud-masked, 
and geolocated Local Area Coverage (LAC) /High 
Resolution Picture Transmission (HRPT) data. The 
Office of Naval Research has supported the development 
of such a data set for most of the Arctic Ocean covering 
the year 1989, selected for best compromise between 
coverage and cloudlessness. Subsequent data are now 
being processed, and will be available through the 
National Snow and Ice Data Center in Boulder, Color- 
ado. Prototyping of data for NASA s planned MODIS 
sensor requires a high resolution AVHRR data set, what 
polar data should be included? The re-analysis and 
archiving of older global data sets are being supported by 
NASA under the Pathfinder program; no high resolution 
polar AVHRR data are presently included. Before 
arguments can be made for better archived polar data 
sets, the utility of such data must be assessed, and a 
justification made. There is something of a chicken-and- 
egg problem here: better data sets must be justified by 
both the scientific needs and maturity of algorithms, yet 
the demonstration of the data’s utility requires improved 
data access. 

To begin to resolve this stalemate, we have initiated a 
review of the status of polar surface measurements from 
AVHRR. As part of this process we have posed the 
following questions, solid answers to which will define 
limitations, strengths, and uses for these data. 

1 . What are the most important polar surface 
parameters that can be derived from the AVHRR 
data, and by what algorithms are they estimated? 

2. How accurately can each be derived? What steps 
could be taken to improve the accuracy? 

3. What validation is needed before these parameters 
can be estimated accurately and routinely for 
scientific research? 

4. What are appropriate sampling objectives? What 
are the optimal spatial resolutions for each of the 
parameters? For what applications are lower 
resolution Global Area Coverage (GAC) data 
useful? For which are higher resolution LAC or 
HRPT necessary? What temporal resolution and 
regional coverage are required? 

5. What difficulties are encountered by users, and how 
can they be mitigated? Are calibration and 
geolocation data adequate? How can the barrier 
of cloudiness be dealt with? Are needed data 
available through current archives? 

A consideration of the first three of these questions 
seemed compatible with this International Glaciological 
Society symposium, so a half-day workshop was sched- 
uled during the meeting. The remaining questions seemed 
more appropriate for a different group of people than 
would be gathered at Boulder. Because of limited time, 
the agenda was focused to the most relevant materials, 
but even then time for discussion was not adequate. For 
that, we apologize to the workshop participants, who 
sacrificed a mountain tour and hike to attend. The 
workshop format consisted of a review of each of seven 


Steffen and others: Snow and ice application of AVHRR 


v 


topics: cloud masking (excluding determination of cloud 
properties), surface temperature, albedo, ice concentra- 
tion, lead statistics, ice motion, and ice sheet properties. 
These seven reviews are presented in the next seven 
sections of this paper. Rapporteurs (the authors of this 
report) reported the conclusions of relevant papers given 
during the regular symposium sessions. 

This report summarizes ideas from the general 
workshop and symposium discussions. To keep the paper 
to a modest length, we have not expanded it into a full 
state-of-the-art review of these seven topics, but such a 
thorough review would be timely. In the final section, we 
draw attention to the most pressing needs for and most 
promising developments of quantitative polar applica- 
tions of AVHRR data. 


2. CLOUD MASKING 

2.1. Introduction 

Cloud masking is of critical importance in polar regions. 
Clouds and snow/ice surface features have similar 
characteristics in both the visible and infrared spectral 
ranges, thereby making polar scene identification ex- 
tremely difficult. Often both high level cirrus clouds and 
low level fog and stratus clouds are optically thin. In such 
cases, retrievals of surface parameters such as albedo and 
temperature are seriously compromised. 

2.2. Algorithms 

As part of the International Satellite Cloud Climatology 
Project (ISCCP), global cloud cover is being retrieved. 
Typical cloud masking algorithms assume that clouds can 
be detected using visible and infrared channel thresholds. 
Reflectance thresholds typically are set about 3% above 
the background, and thermal thresholds typically are set 
about 3°C below the background. Other approaches rely 
upon bispectral thresholding (Minnis and Harrison, 
1984) and a variety of statistical methods (e.g. Saunders 
and Kriebel, 1988). However, Rossow and others (1989), 
Stowe and others (1989) and many others have reported 
difficulties associated with polar cloud cover retrievals. 
Indeed, Landsat imagery shows that clouds often are 
darker than the background snow and ice (Welch and 
others, 1992a). In particular, cloud cover often is 
confused with melt ponds, thin ice and pancake ice. In the 
infrared spectrum, low surface temperatures, strong 
inversions and isothermal structure make cloud discrimi- 
nation difficult. In recent years, it has been shown that 
combinations of spectral and textural measures can be 
effective for polar cloud classification (Ebert, 1987; Key, 
1990a; Welch and others, 1992a). Spectral measures 
typically are based upon AVHRR channels 1, 3 and 4. 
Textural measures typically are computed from channels 
1 and 4; the 3.7 /zm channel (channel 3) is often too noisy 
for reliable textural measures. Difference channels, such 
as chl-ch2 and ch3-ch4, and ratio channels, such as ch3/ 
chi, also can be effective for polar scene identification. 

During the IGS symposium, thirteen groups presented 
a variety of approaches to cloud masking, with various 
levels of success. These techniques can be loosely 


categorized as: (1) spectral thresholding; (2) statistical 
and clustering; (3) textural; and (4) spectral and textural. 
On one extreme, Lindsay and Rothrock (1993) relied 
upon subjective cloud masking. Serreze and others (1993) 
reported problems with the ISCCP cloud detection 
algorithm, and Emery and Fowler (unpublished) tried a 
number of cloud masking schemes, but reported that none 
was definitive. Liger and Fily (1993) used the 1 1 /zm 
thermal and 6.7/zm water vapor channels, but found that 
this combination was not very successful. Maslanik and 
Key (1993) combined AVHRR and microwave SSM/I 
data to improve their retrievals. 

Eppler and Full (1992) applied a Sobel edge gradient 
measure based upon 3 x 3, 5 x 5 and 7x7 convol- 
ution masks, while Comiso (unpublished) used a 
statistical filtering approach. Massom and Comiso (in 
press) used statistical measures based upon ch4 vs ch3-ch4 
for cloud/ice discrimination, ch4 vs. ch2-ch4 for ice type 
discrimination. Aizen and Loktionova (unpublished) 
used channels 3 and 4 not only to discriminate cloud 
cover, but also to retrieve cloud optical thickness and 
aerosol properties. 

Gallegos and others (1993) used cluster shade (3rd 
moment) and data degradation for AVHRR cloud 
masking, and Rees and Lin (1993) applied textures 
generated from the co-occurrence matrix and fractal 
measures to classify Landsat MSS scenes. Welch and 
others (unpublished) applied more than 200 spectral and 
textural measures within an artificial intelligence classi- 
fier. Accuracies greater than 85% were reported with this 
approach. Finally, Chou and others (1993) reported a 
segmentation approach to discriminate cloud and surface 
features. 

Most of the techniques reported during the IGS 
symposium performed adequately on limited datasets. 
However, only one group provided objective analysis of 
classification accuracy. Clearly, before an operational 
polar cloud masking scheme can be validated and 
accepted by the polar research community, it is essential 
that objective measures of accuracy be provided. 
Typically, the Monte Carlo (e.g., Ebert, 1987) and 
Bootstrap (Welch and others, unpublished) approaches 
are used for objective analysis. 

There are three main factors that must be addressed in 
the development of an operational polar cloud masking 
algorithm: (1) the choice of the feature vector, (2) the 
choice of the classifier, and (3) proper identification and 
labeling of regions. First, from the wide range of 
techniques employed by different groups, it can be seen 
that there is no consensus as to what types of signatures 
are needed in a robust polar cloud/surface feature 
classification scheme. However, it is generally recognized 
that spectral information alone is inadequate. Difference 
and ratio channels and textural measures need to be 
explored further to determine if there is an optimum 
feature vector for polar scene identification. In this 
regard, Gray Level Difference Vector textural measures 
seem to offer a good combination of discriminating ability 
with low storage and CPU requirements. 

There is no consensus as to the size of the region over 
which textures should be computed. At present 16 x 16 
pixel regions seem to be preferred. Selection of regions 
that are too small leads to unstable textural measures. 
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However, selection of regions which are too large leads to 
loss of information; i.e., the region is unlikely to contain a 
single class. The 16 x 16 pixel region seems to provide 
the best compromise between textural stability and 
discrimination of features. As a caveat it should be noted 
that textural measures based upon third order (cluster 
shade) and fourth order (cluster prominence) statistics 
generally do not attain textural stability with 16 x 16 
pixel templates; considerably larger template sizes are 
required. Therefore, these measures should be applied 
only after stability analysis is made. 

There may not be an optimum feature vector for polar 
scene identification. Indeed, there seem to be a number of 
combinations of features that are effective for cloud 
masking and classification, a few guidelines are emerging. 
It should be noted that AVHRR channel 3 is particularly 
valuable in discriminating between water and ice surface 
features and between water and ice clouds. Also, 
reflectance values should be divided by the cosine of the 
solar zenith angle. This procedure normalizes reflectances 
from different parts of one scene and from different scenes. 
However, this approach may not be adequate for an 
operational polar scene classification algorithm, Even 
normalized ocean and ice reflectance values are observed 
to vary substantially across a single scan line. This opens 
the question as to whether bi-directional reflectance 
functions must be used to normalize the spectral 
signatures. If so, then a preliminary classification of each 
region first must be made, followed by the application of 
the appropriate bi-directional reflectance function. 
Finally, the solution is iterated to stability. Obviously, 
for large regions composed of mixed classes, this 
procedure can be tedious. Furthermore, the bi-direc- 
tional reflectance functions are known only approxi- 
mately for most classes. 

It should be noted that the cloud classification and 
cloud masking algorithms presented at this symposium do 
not cover the entire range of ideas that are found in the 
literature. For example, it may be possible to utilize 
TIROS Operational Vertical Sounder (TOVS) or High 
Resolution Infrared Sounder (HIRS) data effectively, 
especially for high cloud cover. Another idea which has 
strong merit is the coupling of AVHRR and passive 
microwave retrievals. In spite of their large footprints, the 
microwave channels can provide positive identification of 
surface features. Another approach which has been 
suggested is based upon temporal analysis. It relies upon 
image subtraction to identify regions which have changed 
from one image to another (i.e., cloud motion). Of course, 
one must be careful not to confuse sea-ice motion and 
recent snowfall as cloud cover. Note also that virtually all 
of the methods presented both at this symposium and in 
the literature are effective only during the daylight hours. 
To date an effective “night-time” algorithm has not been 
demonstrated. 

The many groups working on the cloud masking 
problem have used a spectrum of classification schemes. It 
is typical to apply some variant of the Maximum 
Likelihood or Discriminate Analysis approaches. How- 
ever, it was demonstrated at the symposium that artificial 
intelligence (AI) classifiers can significantly increase 
classification accuracy. The main drawbacks of the AI 
schemes are that these methods are not well understood 


by most groups and that the training time may be longer 
than for traditional approaches. On the other hand, in 
operational mode the AI approaches are very fast. It is 
also worth noting that the AI approaches are both 
nonlinear and nonparametric. The advantage is that no 
particular class structure (i.e., normal distribution) is 
assumed; rather, the classifier learns through presentation 
of examples. The AI approaches also require far fewer 
training samples than do traditional schemes. The new AI 
approaches presented at this symposium provide very 
high classification accuracies. The most impressive of 
these techniques is the “Don’t Care” Back Propagation 
neural network (Welch and others, unpublished). Indeed, 
for regions containing pure classes, the accuracy exceeds 
95%. Nevertheless, much. more work is needed to extend 
these results to mixed classes. 

A final issue concerns scene identification and 
labeling. No classification scheme can be expected to 
produce accurate results if the labeling is incorrect. The 
old adage “garbage in, garbage out” is especially 
appropriate here. The analyst needs to examine a wide 
variety of information before labeling a region. A new 
Interactive Visual Image Classification System (IVICS) 
was introduced at the workshop which provides a wide 
variety of analysis tools to the user (Welch and others, 
unpublished). Only using this system have the very high 
classification accuracies (>95%) been attained. The 
analyst needs to take every precaution in labeling; it is 
very easy to make labeling errors in polar scene analysis. 

2.3. Recommendations 

In terms of climate change, cloud classification in terms of 
cloud types, such as stratocumulus, cumulus, and cirrus, is 
somewhat artificial. Future schemes should focus upon 
the retrieval of cloud optical depth, phase (water or ice), 
effective particle size and cloud cover rather than upon 
cloud type. These are the parameters which directly 
impact climate change. Similarly, aerosol size distrib- 
utions and optical depths need to be recovered, at least 
over water surfaces. Knowledge of the optical properties 
of aerosol, diamond dust, and thin water and ice cloud, 
coupled with radiative transfer models, will allow the 
polar research community to retrieve both surface albedo 
and temperature, as well as other parameters, with far 
greater confidence. 

It is recommended that priority be given to the 
development of an operational AVHRR cloud masking/ 
classification scheme during the next three years. 
Additional research needs to be made on: 

1. Computationally efficient and accurate feature 
vectors need to be determined for an operational 
polar cloud masking algorithm. 

2. Coupling of AVHRR, TOVS, HIRS and micro- 
wave (e.g. SSM/I) data to provide positive 
identification of surface features. 

3. Efficient and accurate classification schemes, in- 
cluding AI approaches, have to be tested for pure 
classes, and additional work is needed to extend 
these results to mixed classes. 

4. Bi-directional reflectance functions in polar regions 
(snow and ice surfaces) should be determined and 


applied to normalize the spectral signatures. 

5. Improved visualization and database tools will be 
essential for basin-wide applications. 

As they are developed, these schemes must be 
validated. It is recommended that closer ties be 
established between researchers working on this pro- 
blem; as a start, all interested parties should apply their 
algorithms to a common dataset. Researchers then should 
communicate their findings every six to 12 months to the 
science community to intercompare their results. If 
approached this way, and if validating data from field 
observational programs is made available, we are 
confident that a reliable (>90%) operational polar- 
community cloud masking scheme can be developed in 
three years for daytime imagery. It is too early to estimate 
what accuracy night-time cloud masking approaches may 
achieve. 


3. ICE SURFACE TEMPERATURE 
3.1. Application 

One parameter needed for the estimation of radiative and 
turbulent heat fluxes from satellite data is surface 
temperature. Little effort has been directed to the 
retrieval of ice surface temperature (1ST) in the Arctic. 
The reason is not one of methodology, but rather our 
limited knowledge of atmospheric temperature, humidity, 
and aerosol vertical, spatial, and temporal distribution, 
the micro-physical properties of polar clouds, and the 
spectral characteristics of snow, ice, and water surfaces. 
Over the open ocean the surface is warm, dark, and 
relatively homogeneous. This makes sea surface tempera- 
ture (SST) retrieval, including cloud clearing, a fairly 
straightforward task. Over the ice, however, the surface 
within a single AVHRR pixel is likely to be highly 
heterogeneous, with the spatially integrated emissivity 
dependent upon the exact mixture of ice of various 
thicknesses, open water, and snow cover. In addition, 
Arctic cloud amounts are high (typically from 60-90%) 
complicating 1ST retrieval even further. 


3.2. Algorithms 

The general approach to estimating surface temperatures 
from thermal channels is to relate satellite observations to 
surface temperature observations with a regression model. 
Massom and Comiso (in press) have followed this 
approach regressing brightness temperatures from the 
AVHRR against station temperatures in the Bering and 
Greenland seas. The estimated surface temperatures have 
an accuracy of 1-3 K. 

Lacking a sufficient quantity of corresponding satellite 
and surface observations, however, satellite radiances or 
brightness temperatures can be modeled by forward 
calculation of the radiative transfer equation. This 
approach is commonly used for SST retrieval (e.g., 
McClain and others, 1985). SST retrievals with two 
"split-window” channels have been shown to have an 
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absolute accuracy as high as 0.5-1 K (RMS error). Land 
surface temperature (LST) estimates commonly follow 
the same approach but are generally less accurate due to 
the larger variability of surface conditions, and errors of 
2-3 K are common. This method of SST retrieval has 
been employed for 1ST retrieval from AVHRR data by 
Key and Haefliger (1992). Temperature and humidity 
data from Soviet ice islands, skin temperatures estimated 
with an energy balance model, and modeled directional 
emissivities of snow were used in the simulation of 
AVHRR radiances. A regression model was developed 
that includes the two split-window channels and satellite 
scan angle. RMS error in the estimated ice surface 
temperature is stated to be less than 0. 1 K in all seasons. 
However, this value is based on the regression analysis of 
the modeled data only; it does not include measurement 
error, emissivity differences due to mixed pixels, or 
unmodeled atmospheric effects such as ice crystal 
precipitation (ICP). In fact, Lindsay and Rothrock 
(1993), who have applied the same coefficients to Arctic 
AVHRR data, report the AVHRR ISTs to be biased by 
about + 4K in January 1989 based on comparisons with 
drifting buoys. They speculate that this large an error can 
be attributed primarily to the presence of "diamond 
dust” (ICP). 

The basic SST approach described above was also 
taken by Haefliger and others (1993) with AVHRR data 
over the Greenland ice sheet. Temperature and humidity 
profiles from a research site on the Greenland ice sheet 
were used in the analysis. Skin temperatures were 
estimated based on upwelling longwave radiation. These 
same data were used for validation. Surface temperatures 
were estimated to within 0.5 K in the worst case, with a 
mean of about -0.2 K. 

While not utilizing AVHRR data directly, Maslanik 
and Key (1993) presented preliminary results on the use 
of SSM/I passive microwave data for 1ST estimation. The 
coefficients from Key and Haefliger (1992) were used to 
estimate the surface temperature within 25 km cells in 
winter and then converted to ice temperatures assuming 
SSM/I-derived sea-ice concentrations. These ice tem- 
peratures were then used to derive an "effective 
emissivity” for each SSM/I channel, which relates the 
physical temperature at some frequency-dependent depth 
within the snow/ice pack to the skin temperature. While 
errors up to 4K were reported, greater accuracy would 
probably result if the effective emissivities were deter- 
mined for smaller regions with more homogeneous 
emissivities and if SSM/I swath data were used rather 
than the daily averages. 

Finally, although not discussed during the conference, 
the cloud product of the ISCCP contains global surface 
temperature estimates. For the polar regions these are 
based on AVHRR data. At the time of this writing, data 
are available from mid- 1983 through 1990. The proced- 
ure involves a radiative analysis for each pixel of a 
reduced-resolution product, and therefore differs from the 
regression methods described above. Temperatures 
appear to have a positive bias of a few degrees in 
winter, although this will probably improve when the 
original data are reprocessed with a special polar cloud 
detection algorithm (personal communication from W. 
Rossow, 1992). 
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3.3. Recommendations 

The difficulties in retrieving 1ST from AVHRR thermal 
sensors are not related to methodology per se, but rather 
to our inadequate knowledge of the atmospheric and 
surface conditions encountered in the polar regions. Are 
current estimates accurate enough calculating energy 
fluxes over ice-covered waters for monitoring climate 
change? In terms of modeled accuracy only, 1ST can be 
retrieved to within 0.1— 0.2 K. A more realistic estimate, 
which includes uncertainties in the measurement and 
physical (atmospheric and surface) systems, is probably 
1-4 K. This translates into an error in upwelling 
longwave flux of 3-1 1 Wm 2 at 230 K or 5—18 Wm at 
270 K. These values are of the same magnitude as 
expected changes that we need to detect, so that 
improvements are certainly required. 

In order to obtain more accurate estimates of 1ST 
over the polar regions we therefore suggest the following: 

1. Cloud masking methods (including detection of ice 
crystal precipitation) should be improved and 
automated over snow and ice surfaces, especially 
during the dark months. 

2. Additional information and better understanding of 
the spatial and temporal distributions of ice crystal 
precipitation and aerosols are needed. This parti- 
cular atmospheric phenomenon presents perhaps 
the greatest uncertainty in 1ST retrievals from 
thermal data during late fall, winter, and early 
spring. It is difficult to detect by its spectral 
characteristics alone. 

3. During summer the water emissivity value should 
be used mainly to reflect the radiative properties of 
melting ice, melt ponds and open water. During 
winter 1ST retrieval errors due to the emissivity 
difference of snow and water are less than 0.2 K and 
therefore can be neglected at a first approximation. 

4. While a considerable amount of surface tempera- 
ture data has been collected by aircraft during 
various field programs, experiments specifically 
designed for the validation of satellite remote 
sensing algorithms are needed to cover the range 
of surface and atmospheric conditions found in the 
polar regions. 

4. NARROW-BAND ALBEDO 
4.1. Introduction 

Satellite-derived albedo values are useful for the estima- 
tion of the shortwave radiation balance, an important 
term in surface energy balance studies. For the study of 
spatial and temporal changes of surface albedo, satellites 
provide the only data source due to the inherent dynamic 
change in the sea ice cover. Over ice sheets, surface 
reflectance is a good indicator for classifying different 
surface types such as glacier ice and various snow facies. 
Snow and ice surfaces have relatively high albedo values 
(50-90%), whereas the water surface has one of the 
lowest albedo values (5-10% for high solar elevations). 
Snow and ice albedo depend on the solar incidence angle, 


snow grain size and geometry, impurities in the snow, 
water content in the snow cover, surface topography, and 
reflection angle. The albedo should be distinguished from 
the reflectivity, which refers to one specific wavelength. 

4.2. Algorithm 

The AVHRR multispectral narrow-band radiometric 
scanners with small fields of view need proper adjustments 
so that their data can be used to estimate top-of-the- 
atmosphere albedo values (planetary albedo). The 
conversion of shortwave narrow-band data into average 
planetary albedo requires adjustment of the spectral 
expansion from filtered narrow-band to unfiltered broad- 
band radiance of the shortwave spectrum from 0.24 to 

4.2 fim, The planetary albedo for snow surfaces using 
AVHRR narrow-band satellite data can be derived 
through a multivariant regression analysis of Nimbus-7 
earth radiation broad-band data and AVHRR narrow- 
band data (Wydick and others, 1987). This statistical 
approach was done for NO A A- 7, and still needs to be 
determined for the NOAA satellites 8 through 12. 

To derive surface albedo values from AVHRR 
narrow-band radiance measurements, the scattering and 
absorption of the atmosphere has to be included. This was 
successfully demonstrated by Koepke (1989) through 
numerical simulation that used a linear relationship 
between clear sky planetary albedo and surface albedo. 
The coefficients he derived are given as a function of solar 
zenith angle, optical depth of the atmosphere due to 
aerosols, integrated ozone content of the atmosphere for 
AVHRR channel 1, and water vapor content for 
AVHRR channel 2. An alternate way to derive surface 
albedo values from AVHRR radiance measurements is to 
use a radiative transfer model such as LOWTRAN to 
correct for scattering and absorption as demonstrated by 
Haefliger and others (1993). 

4.3 Modeling and experimental measurements 

The knowledge of temporal changes in the spectral signal 
of snow is important for the understanding of albedo 
changes. The spectral albedo of snow in the solar 
spectrum is strongly a function of both the optically- 
equivalent grain size, and concentration of absorbing 
impurities such as soot. Davis and others (1993) have 
shown, in their statistical analysis derived from stereo- 
logical measurements, that the optically-equivalent grain 
size is a sphere with equal surface- to- volume ratio as the 
ice matrix. They used a one-dimensional energy and mass 
budget model that interacts with the snow cover and used 
outputs from this model as boundary conditions for a 
discrete ordinate radiative transfer model to predict the 
spectral bidirectional reflectance of the snowpack. These 
modeled results have still to be verified by in situ 
measurements. 

Bidirectional reflectance of experimental measure- 
ments were compared with model results from a 48- 
streams discrete-ordinate model (Steffen and Schweiger, 
unpublished). The model follows the Wiscombe and 
Warren (1980) approach where the two-stream approach 
is replaced with a discrete-ordinate solution (48 streams) 
to the radiative equation in order to allow modeling of the 
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Bidirectional Distribution Reflectance Function (BDRF). 
The model assumes Mie-scattering of optically equivalent 
ice spheres. Small amounts of carbon soot had to be 
included in the medium to achieve a match between 
measured and modeled reflectances in the visible region. 
Comparison of model results with field measurements 
showed, in general, a good agreement if a soot amount of 

0.30 parts per million weight (ppmw) was used. If the 
model run was performed with 0.00 ppmw soot amount, 
the spectral reflectance in the visible range (0.3-0. 8 /nn) 
was approximately 0.1 higher than the in situ measure- 
ment. 

Experimental measurements of optical properties were 
carried out in the laboratory for 70 snow samples 
collected in the Alps (Sergent and others, 1993). The 
carbon soot content was measured for all snow samples. 
The spectral reflectance between 0.4— 1.0/im was ob- 
tained with a monochromator and comparison with the 
Wiscombe and Warren (1980) model showed good 
agreement. The study showed that crystal types and 
sizes (usually described in a subjective and simplified 
manner) have to be determined with great accuracy (i.e. 
statistical analyses of macro-photographs) to show a 
useful relation between snow grain-size, impurity content, 
and spectral reflectance. 

Surface shortwave radiative fluxes over polar oceans 
can be obtained on a climatological basis from joint use of 
AVHRR and SSM/I data (Kergomard and others, 
1993). The directional AVHRR radiances at the top of 
the atmosphere were simulated with a radiative transfer 
model with realistic atmospheric conditions, and com- 
pared in an iterative mode with average AVHRR visible 
radiances for cloudy pixels to derive optical thicknesses of 
the cloud cover. The SSM/I data provided the sea-ice 
concentrations and consequently, the surface reflectance. 
From the retrieved optical thickness, the incoming, 
absorbed and reflected solar fluxes (0.2-2. 8 pm) were 
computed by using a two-stream narrow-band model. 
The results agreed well with in situ measurements, but 
they are strongly dependent on cloud base temperature. 

4.4. Errors 

Prelaunch radiometric sensitivity of AVHRR visible and 
near infrared channels (1, 2) are determined with an 
integrating sphere by the manufacturer of the sensors to 
indicate their precision (Abel, 1990). The integrating 
sphere, however, has also a degrading sensitivity with 
time. Therefore, prelaunch calibrations, which are based 
on the intercomparison with the integrating sphere, will 
always have a relative error of up to a few per cent (i.e. 
the sphere calibration showed a drop in radiance of 10- 
16% between 1983 to 1987 (Abel, 1990)). 

A second and probably more important error source is 
the rate of degradation of the sensors over time. For the 
visible channel of NOAA 9, degradation rates of up to 6% 
per year were found (Staylor, 1990). Haefliger and others 
(1993) derived degradation rates of 0.6% per year for 
channel 1 and 4% per year for channel 2 for NOAA 1 
based on a single comparison of in situ measured and 
modeled (AVHRR) reflectance values from the Green- 
land ice sheet. 

In situ measurements of water vapor, aerosols, and 
ozone distribution with height are usually not available 


and must be parameterized. This parameterization bears 
another source of error which is difficult to estimate. 
Aerosol content can vary with time and space signif- 
icantly in Arctic regions which can reduce the reflectance 
by up to 10% as reported by Lindsay and Rothrock 
(1993). 

Satellite sensors with narrow field-of-view scanners 
measure the bidirectional reflectance at one or a few 
angles. If the reflecting surface, in our case snow, ice and 
water, is not Lambertian, than the measurements have to 
be corrected by the anisotropic reflectance factor (ARF), 
which is defined as the bidirectional reflectance normal- 
ized relative to the hemispheric reflectance. The first 
airborne bidirectional radiances of snow-covered surfaces 
have been reported by Hall and others (1993). For wet 
snow surfaces, the ARF at large solar zenith angles can be 
as large as 2 (Steffen, 1987). For new snow the ARF is 
close to 1, and most likely can be neglected for AVHRR 
applications. 

4.5. Recommendations 

The overall accuracy of AVHRR-derived albedo values 
depends on a number of error sources, such as: pre-launch 
calibration, sensor degradation, uncertainty in aerosol, 
ozone and humidity concentrations over time and space, 
anisotropic reflectance of snow and ice, cloud contamina- 
tion of pixels, and conversion from narrow-band albedo 
to broad-band albedo values. With the given uncertain- 
ties the expected error in AVHRR-retrieved surface 
albedo could be as large as 10%. In a simple energy 
balance model this error (A albedo = + 10%) translates 
to a 3°C reduction in air temperature on an annual mean. 
This uncertainty is too large to be useful in climate 
sensitivity studies, and therefore, the following improve- 
ments are recommended. 

1. AVHRR radiometric degradation has to be 
monitored carefully over time (no onboard calibra- 
tion). Large scale reference areas of well known 
spectral reflectance properties, such as deserts or 
dry snow areas of ice sheets, could be used for 
frequent calibrations. 

2. Conversion of planetary albedo to surface albedo is 
very sensitive to humidity, aerosol, and ozone 
concentrations in the atmosphere. In situ measure- 
ments of these variables are needed to characterize 
the atmospheric conditions to correct for absorption 
and scattering (two-stream models). 

3. A statistical climatology of aerosol, ozone, and 
humidity distribution derived from satellite soun- 
ders is needed for different Arctic regions to reduce 
the uncertainty of large scale albedo retrieval. 

4. Radiative transfer models with high spectral 
resolution should be used to correct for absorption 
and scattering based on statistical climatology as 
mentioned above (3). 

5. ICE CONCENTRATION 
5.1. Introduction 

For some time AVHRR data have played an integral role 
in general ice-forecastinp’ oneratinns that nrnHnrp mam nf 
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the ice edge and chart the location of large polynyas and 
wide leads (Hufford, 1981). Although effective automated 
means have not been developed to extract ice concent- 
ration from the AVHRR data stream, analysts who 
prepare operational ice charts through conventional 
manual interpretation methods use AVHRR images to 
map the ice edge on a routine basis. 

More extensive use of AVHRR data to map ice 
concentration within the pack and to extract ice type 
information, whether through manual analysis, subjective 
classification techniques, or objective automated routines, 
has been limited, in large part by two factors. First, 
discrimination of thin ice from open water is next to 
impossible, particularly in the vicinity of open leads and 
polynyas where significant atmospheric water vapor and 
ice crystals bias radiances sensed from space. Second, 
interpretation of pixels that include mixtures of ice types 
that occur at subresolution scales, such as in leads or in 
areas where brash or small floes are common, is 
ambiguous. 

5.2. Algorithm and accuracy 

AVHRR ice edge locations for the Bering Sea in March 
1988 were validated with Landsat imagery by Emery and 
others (1991a). The AVHRR ice edge location was found 
using a subjective threshold method to classify between 
water and ice. Results show the AVHRR edge to be 
consistently seaward of the Landsat edge. Mean deviation 
is 1.8 km in channel 2 data, and 2.9 km in channel 4 data 
with rms error of 3.2 km and 4.4 km. 

Derivation of ice concentration from AVHRR 
imagery requires that image pixels be classified as either 
water, ice, or cloud. Continuing progress in development 
of cloud masking algorithms as discussed in Section 2 will 
effectively reduce the problem to discrimination between 
ice and water in the near future. For clear cloud scenes, 
ice concentration can be achieved during winter months 
when good thermal contrast exists between ice and water 
surfaces. 

Issues related to the derivation of ice concentration 
and ice type information from AVHRR imagery were 
addressed by two papers presented during this sym- 
posium. The first paper, by Key and others (1993), 
established a basis for evaluating the likelihood of 
observing a particular feature if its size with respect to 
the sensor beam spot, its contrast with adjacent ice, and if 
the conditions of the atmosphere are known. The authors 
base their analysis of contrast and atmospheric conditions 
on numerical simulations that use the LOWTRAN-7 
model. Empirical approaches that degrade fine resolution 
Landsat images of leads to the coarser AVHRR 
resolution form the basis for analysis of scale relation- 
ships. Whereas work reported by Key and his co-authors 
was developed specifically to assess potential error in lead 
statistics derived from AVHRR imagery, the results have 
broader implications with respect to feature detection in 
general and with respect to derivation of ice type and ice 
concentration information. 

The second paper, presented by Massom and Corniso 
(in press) demonstrated that frequency distributions 
formed by values of pixels in an AVHRR image are 
nnlvrrmHal and carrv information about the aerial extent 


of different ice types that are present in the scene. Massom 
and Corniso assign ice types to individual modes by 
interpreting relative brightnesses and temperatures that 
define pixels within the modes. The areal extent of each 
ice type in the scene then can be estimated as a function of 
the number of pixels within each mode. The method 
depends, of course, on reliable cloud-masking procedures 
to identify cloud pixels and remove them from consider- 
ation. 

These two studies complement each other and, if 
future efforts are coordinated, could form the basis for 
deriving improved ice concentration estimates from 
AVHRR data. Massom and Corniso (in press), on one 
hand offer the potential for ice classification on a pixel-by- 
pixel basis. The weakness of this approach is two-fold. 
First, pixels that include small features at subresolution 
scales (e.g., leads narrower than 1 km in the case of 
AVHRR nadir pixels) will be misclassified because their 
brightness represents the integrated radiance of more 
than one ice type — ice within the feature and ice 
adjacent to the feature. Second, variability of the 
atmosphere with time and location changes radiances 
measured from space and shifts the position of modes 
within the frequency distribution. Thresholds that define 
boundaries between ice types thus vary with respect to 
location within a single scene. Key and others (1993), on 
the other hand, provide the means to address these issues, 
to assess their impact on error inherent in AVHRR ice 
concentration products, and to adjust derived concentra- 
tion estimates by appropriate values so that they reflect 
actual surface conditions more accurately. For example, 
frequency distributions used by the Massom-Comiso 
method can be modeled to determine the extent to 
which ice concentration retrievals are robust with respect 
to specific atmospheric profiles, varied assemblages of ice 
types, and collections of subresolution features (leads, 
floes, brash ice, and bodies of slush and shuga) that occur 
at different scales. Lead pixels, which can be identified 
using segmentation techniques and Hough transform 
methods (Duda and Hart, 1972; Fetterer and Holyer, 
1989) and other pixels that are likely to include multiple 
surface types can be analyzed independently so as not to 
increase overall error. These mixed-surface pixels, which 
are most likely to include thin ice, can be analyzed further 
to narrow the range of possible surface types and feature 
widths required to produce the observed brightness 
signature. 

5.3. Recommendations 

Progress in three areas is requisite to utilize fully the 
results of the two studies mentioned above (Massom and 
Corniso, cited above; Key and others, 1993) and to 
improve ice concentration algorithms in general. Pre- 
requisites for AVHRR-derived ice concentrations based 
on threshold methods are accurate ice surface and 
reflectance measurements (see Sections 3 and 4). There- 
fore, programs that provide for full characterization of 
radiometric properties of both sea ice and the atmosphere 
at AVHRR frequencies, improved methods of modeling 
the surface temperature field, and implementation of 
cloud-masking algorithms are needed. 
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1. Ice properties . Optical and infrared signatures of 
different ice types and snow types must be described 
in terms of reflectance (possibly bidirectional 
reflectance), and thermal emissivity in order to 
define radiometric contrast for the work of Key and 
others, (1993), and to verify ice interpretations 
made in the course of the work of Comiso 
(unpublished). Limited data that describe some of 
these parameters have been published, both with 
regard to albedo of some types of sea ice (Grenfell 
and Perovich, 1984; Perovich, 1991) and bidirec- 
tional reflectance of snow (Dirmhirn and Eaton, 
1975; Kuhn, 1985; StefTen, 1987). Additional 
measurements are needed to convert from broad- 
band albedo to narrow-band albedo as measured 
by the AVHRR. Further, information is needed to 
characterize the transition from open water to thin 
ice, which is essential for monitoring heat and brine 
fluxes. 

2. Characteristics of the atmosphere. Data that define 
atmospheric characteristics (water vapor, ice crys- 
tals, diamond dust, temperature) in polar regions, 
and thus the magnitude of bias these parameters 
contribute to space-based optical and infrared 
measurements of surface temperature and bright- 
ness, are not collected routinely except at widely 
scattered shore stations. Temporal and spatial 
variability is thought to be significant but is poorly 
documented at best and, in the central Arctic basin 
and in ice covered Antarctic waters, day-to-day 
variability is virtually unmonitored. Effort should 
be directed toward developing improved methods 
to model scattering, emission, and attenuation 
characteristics of the atmosphere using satellite 
data. Emphasis should be placed on innovative use 
of new sensors that will fly on EOS-era platforms. 
In particular, opportunities to employ space-based 
laser sensors as scatterometers to infer the presence 
of diamond dust should be explored. 

3. Implementation of cloud detection algorithms . Automated 
ice concentration algorithms will require accurate 
cloud masking procedures. Progress toward effec- 
tive cloud classification algorithms has been 
significant over the past five years. Efforts pursued 
by disparate groups of researchers should be 
coordinated to produce, within three years, a 
hybrid algorithm that embodies the best aspects of 
each approach, with transition to an operational 
system within five years. 

6. LEAD STATISTICS 
6.1. Application 

Mapping leads in AVHRR imagery is a task with several 
applications. Heat flux through leads plays an inordi- 
nately large role in the Arctic heat budget. Lead width, 
lead density, the presence or absence of new ice in leads, 
and the orientation of leads relative to the wind are all 
important parameters when calculating surface fluxes. Ice 
rheology is revealed through large- and small-scale lead 
patterns. The opening and closing of leads is an effect of 


ice dynamics, which in turn affects the physical and 
biological oceanography of the upper ocean. Finally, the 
operational community would like to improve map and 
forecast lead width, orientation, and spacing, as these 
parameters influence the navigation of ships and 
submarines. 

6*2. Analysis 

The advantage of AVHRR for lead mapping is its 
synoptic view of Arctic basins and its frequency of 
coverage, which allows the study of lead formation in 
response to synoptic weather patterns and of regional and 
seasonal variations in lead distributions. Mapping leads is 
complicated by the same factors which challenge the use 
of AVHRR for ice concentration or surface temperature 
(see Sections 1-5). An added difficulty is that lead widths 
are generally assumed to follow a power law distribution. 
Relatively few leads are wide enough to fill the sensor field 
of view. The sensor resolution is about 1 km at nadir, and 
degrades to about 4 km at the scan limbs. Therefore 
virtually all “lead” pixels are mixtures of lead and 
surrounding ice, and the majority of leads are too narrow 
to be detected. 

Walter and Overland (1993) illustrate that not all 
applications need wait for the full resolution of the issues 
above. Using AVHRR imagery from March and April 
1991, they relate changes in lead patterns to wind fields. 
A long lead running diagonally across the Beaufort Sea 
revealed by sequential imagery acted as a slip boundary. 
North of this boundary, a diamond-shaped pattern of 
leads was present on two scales. At both scales, the angle 
of intersection between leads was about 30°. This pattern 
supports the theory of Erlingsson (1988), whose analysis 
of the mechanism of ice failure points to self-similarity of 
lead patterns. South of the slip boundary, a diamond- 
shaped pattern is superimposed on a previously existing 
arch-shaped pattern. The dual pattern is related to 
changes in the wind field relative to the coastline and slip 
boundary to the north. This work emphasizes the 
feasibility of deducing ice mechanics from observable 
lead patterns, and the importance of the location relative 
to the direction of forcing of the slip, no-slip, or shear 
boundaries which confine lead pattern regimes. 

Whether imagery is interpreted manually, as in the 
above studies, or automatically, some determination must 
be made of the conditions under which leads can be 
detected, with a view towards error analysis. In thermal 
data, the detectability of a lead is a function of the 
temperature contrast between a lead pixel and the 
surrounding ice, the size of the feature relative to the 
field of view of the sensor, and the effect of the intervening 
atmosphere. Key and others (1993) express detectability 
as a contrast ratio of the brightness temperature of leads 
measured at the sensor to that of thick ice. The ratio is a 
function of measurement scale (the ratio of the width of 
the lead to the sensor field of view) and of atmospheric 
optical depth. To model the function, the effect of 
measurement scale is obtained through degrading Land- 
sat imagery of leads, while the effect of optical depth is 
obtained through modeling surface emissivities and 
radiative transfer using historical data. One result of the 
investigation into measurement scale is the observation 
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that wide leads “grow” as imagery is degraded, while 
narrow leads disappear as averaging quickly reduces their 
contrast with surrounding ice. This implies that lead- 
width statistics from AVHRR imagery should be adjusted 
for this effect. In investigating the effect of optical depth, 
the authors find that while brightness temperatures 
decrease slowly for four modeled ice thicknesses as 
aerosol optical depth increases, the brightness tempera- 
tures quickly converge with rising ice crystal optical 
depth. When ice crystal optical depth is greater than 3.0, 
it is difficult to distinguish open water in leads from ice. 

Key and his coauthors (1993) supply a rational basis 
for selecting a threshold to segment leads in digital 
thermal imagery. An alternative to a straight threshold is 
the potential open water concept of Lindsay and 
Rothrock (1993). In this method, tie points for the 
reflectance or temperature of open water and of thick ice 
are used to establish the potential open water in each 
pixel. By applying a threshold to define lead pixels based 
on potential open water, no assumption need be made 
concerning the ice thickness or mixture of ice and open 
water in a pixel. The authors then derive lead widths 
using randomly placed transects and find that lead (and 
floe) widths are described by a negative exponential 
distribution over a range of regions and seasons. Lead 
width distributions obtained in this manner can serve as 
input to climate models. 

Lindsay and Rothrock manually screened clouds to 
obtain data sets for studying regional distributions. This 
method is accurate but can be time consuming and can 
result in the loss of good data if the image interpreter is 
too conservative. However, constructing an algorithm for 
cloud detection in the Arctic is challenging (see Section 
2). Rather than search for clouds directly, Eppler and 
Full (unpublished) developed a method of automatically 
distinguishing areas covered by leads from other surface 
types by applying edge filters to enhance the spatial 
signature of leads, and then classifying feature vectors 
drawn from cells in the filtered image. Lead-covered cells 
can then be supplied to a lead-statistics algorithm. 

6.3. Error analysis 

Errors in digital lead mapping are errors in detection (or 
how many lead pixels are missed), and errors in 
placement (or in the derived value of orientation, 
width, area, number, and spacing). Once lead pixels 
have been identified, errors in placement depend entirely 
upon the algorithm used to derive statistics. Key and 
others (1993) laid the groundwork for evaluating errors in 
detection when a thermal threshold was used to segment 
lead pixels by answering the question “What is the 
expected contrast between lead and ice for a given lead 
width and atmospheric condition?”. If the variability of 
the temperature of background ice, the distribution of 
lead widths, and the variability of optical depth can be 
estimated, then the possibility of missed or falsely detected 
pixels can be evaluated. The change in field of view across 
the swath, and the effect of mapping from swath to 
gridded data, must also be taken into account. 

6.4. Recommendations 

tv m™r lead maDDine work has relied on manual 


analysis of imagery (e.g. Barry and others, 1989; Marko 
and Thomson, 1977) to yield information on ice rheology 
or lead formation. Using lead statistics as input to models 
has been largely unexplored, and is an application which 
puts more exacting requirements upon lead mapping 
techniques. Appropriate digital data sets and software 
tools can facilitate this work (see, for example, Lindsay 
and Rothrock, 1993), as will progress in the following 
areas: 

1 . Automated algorithms. The success of lead statistics as 
input to models will depend upon accurate 
climatologies. Automated analysis methods are 
needed to handle the volume of data required for 
good climatologies. Although not discussed at the 
Boulder symposium, two publications describe such 
algorithms. Key and others (1990b) use a rule- 
based system to decide which features constitute 
leads and to produce distributions of width, length 
and orientation. One advantage is that rules for 
lead definition can be fitted to the application. 
Fetterer and others (1990) extract distributions of 
lead width, orientation, and spacing using a Hough 
transform-based line finding technique. Rather 
than depending solely on reflectance or thermal 
contrast to identify lead pixels, algorithms can use 
the linearity of leads. This is an important 
advantage given the variability in lead contrast. 
Of course, errors in derived lead parameters due to 
the algorithm itself must be accounted for. This can 
be an onerous task when the algorithm involves 
many processing steps, and was not performed in 
the cited work. 

2. Lead width distributions . Two opportunities for 
extending AVHRR measurements to smaller 
scales were put forth by symposium participants. 
The scale-invariant nature of lead patterns ob- 
served by Walter and Overland (1993) suggests 
that if the boundaries of pattern regimes are 
identified, smaller scale patterns can perhaps be 
inferred. Key and others (1993) attempted to use 
lead width distribution from AVHRR to predict 
the true (full-scale) width distribution. Results were 
overly sensitive to the AVHRR data used as input, 
but accuracy could be improved using coincident 
SAR or Landsat data to fix the predicted 
distribution at the high end. More comparisons of 
AVHRR with higher resolution data are called for 
in order to understand better the relationship 
between the pattern of large and small scale leads 
and the variability in lead width distributions. 

3. Access to more data. One year (1989) of a five-year 
data set of digital AVHRR data covering the entire 
Arctic as frequently as every three days has been 
assembled under the sponsorship of the Office of 
Naval Research (Fetterer and Hawkins, 1993). The 
data set was compiled to address the problem many 
researchers have in finding sufficient quantities of 
processed, cloud-free data over a wide range of 
seasons and regions. Broader use of this data set for 
lead studies and more detailed correlation analysis 
with Landsat, Spot and other high resolution 
images are encouraged. 


Steffen and others: Snow and ice application of AVHRR 


7. SEA-ICE MOTION 

7.1. Applications 

Ice motion measured from sequential satellite imagery is 
applicable to a range of geoscientific problems. In the 
momentum balance for sea ice, the water stress term 
depends directly on ice velocity, and the divergence of 
internal ice stress depends on the spatial gradient of ice 
velocity. The ice mass balance equation has an advective 
term, given by ice velocity. New ice production and the 
resulting brine flux to the ocean surface depends on 
spatial differences of ice velocity or, for shore leads and 
polynyas, on the offshore component of ice velocity. The 
piling of ice mass into ridges depends on convergent and 
shearing ice motions. The ice velocity field is a strong 
boundary condition (certainly more directly observable 
than stress) on the ocean surface. Ice velocity data 
enhance operational forecasts of polynya and lead 
openings or closings, and of the advection of multiyear 
floes into navigation lanes. The motivation for acquiring 
ice motion data is high. 

7.2. Algorithms 

Ice tracking consists first of identifying the same ice 
feature in each of two sequential images, A and B, and 
second, of specifying the geographic location of these 
points in the two images. Each of these two steps has an 
algorithm. The algorithm most widely used for tracking 
ice in AVHRR imagery is maximum cross-correlation. 
This technique finds a local maximum in the cross- 
correlation as a function of the spatial lag (in x and y) 
between the two scenes (Ninnis and others, 1986; Emery 
and others, 1991b; Zhang, 1991). The spadal lag of the 
peak correlation indicates the distance the ice has moved. 
Edge feature tracking and floe boundary tracking 
algorithms applied recently to SAR ice tracking 
(McConnell and others, 1991; Banfield, 1991) have 
generally not been applied to AVHRR. 

Geolocation is accomplished by a second algorithm. 
The procedure is to adjust parameters in an orbit model 
to give correct coordinates of known ground points, and 
to use the new orbit information to compute new 
locations for all points in the image (Baldwin and 
Emery, 1993; Emery and others, 1989). When compar- 
ing two images as when tracking ice, one image is 
sometimes “nudged” (translated uniformly) so that fixed 
(land) portions of the image more perfectly overlay each 
other. The procedure of warping one image to correspond 
to another, which is common for overlaying, say, Landsat 
scenes to maps, is inappropriate for AVHRR sea-ice 
tracking; in sea-ice scenes with land in a portion of the 
scene, warping can cause severe distortions away from 
land. 

Three papers presented at the symposium dealt with 
ice motion or image geolocation. Heacock and others 
(1993) discuss the difficulties of tracking seasonal ice off 
the coast of Newfoundland and Labrador. In addition to 
the usual problem of cloud cover, floe sizes well below the 
resolution of AVHRR cause difficulty to both correlation 
and feature tracking algorithms; the tracker must 
recognize the patchiness in the ice cover. Consequently, 


tracking works best for quite short time intervals, 6 to 48 
hours. Feature identification errors average 2.5 km. Using 
a “nudging” navigation technique, they obtain geoloca- 
tion errors of about 2.7 km. 

Emery and Fowler (1992) reviewed an operational 
version of a correlation ice tracker developed for the 
Navy/NOAA Joint Ice Center (Emery and others, 
1991b). The algorithm adds a spatial filtering techni- 
que, eliminates Fast Fourier Transforms to compute 
correlations, and allows for irregular correlation areas 
near coastlines. Ice velocities computed from different 
AVHRR channels (2 and 4) show almost no difference. 
The use of low-resolution (simulated GAC) imagery 
provides roughly one quarter of the useable motion 
vectors found with the higher resolution LAC images, so 
the field of motion loses much spatial detail. Clearly, 
Fram Strait and the Greenland Sea are ideal for AVHRR 
ice motion observations: the ice motions are large enough 
to swamp the errors, and there is land on both sides of the 
strait to maximize the geolocation accuracy. 

Baldwin and Emery (1993) describe a geolocation 
algorithm that computes from two or more ground 
control points both the satellite orientation (roll, pitch, 
and yaw) and the satellite position at a known time (the 
time the image is acquired). The algorithm can compute 
both the earth locations of image pixels (“satellite 
perspective” or orbital data) and resampled pixels on a 
chosen geographic projection (“indirect navigation” or 
geolocated imagery). The geolocation accuracy of the 
method is 1 km. The paper describes the steps required in 
the calculation, and techniques used to make the 
computing more efficient. 

7.3, AVHRR compared with other motion data 
sources 

Ice motion data are also available from drifting buoys 
and from other satellites including those with synthetic 
aperture radar (SAR). In light of SAR’s substantially 
higher resolution and lower geolocation errors, one can 
ask what role AVHRR should play in a routine ice 
tracking capability. At present AVHRR coverage has 
some advantages to that of SAR. AVHRR has a wide 
swath (about 2750 km); it provides an occasionally 
excellent basin-scale perspective. AVHRR HRPT is 
received at many ground stations giving a great deal of 
Arctic data, but only some 10 to 30% is cloud-free. The 
ERS-1 and JERS-1 SAR satellites provide narrow swath 
coverage (100 km) leaving many areas unseen. Their 
pixel spacing is far superior to AVHRR’s 1 km or greater 
— 12.5 m at high resolution and 100 m at low resolution. 
The Optical Line Scanner on the DMSP satellites with a 
600 m pixel spacing is another excellent tool for testing 
the accuracy of ice motion data. The advangtage of buoys 
is the length and continuity of their data record; their 
positional accuracy (equivalent to <r g ) is several hundred 
meters, sufficient to provide valuable comparisons with 
AVHRR ice motions. 

7.4. Errors 

Errors in the ice motion are described as uncertainty in 
the eeoeraohic location of anv nix^l in pitVipr crAnA* 
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uncertainty in identification of the same feature in both 
scenes. Although error models are generally lacking in the 
literature, both the geolocation error and the feature 
identification error have apparently been taken to be 
zero-mean, Gaussian errors with standard deviations a g 
(for geolocation) and Oi (for feature identification). We 
take discussion of “errors” in the literature to refer to 
these two standard deviations. The geolocation error 
applies to the feature position in each image; the feature 
identification error applies to the feature position in 
image B, the feature being defined by its position in image 
A. Although many geophysical applications require 
velocities, the actual measurement is of position and its 
temporal and spatial differences. A component of 
displacement between images A and B is 

u = (xb - za) (1) 

and has an error with a standard deviation, assuming 
uncorrelated geolocation errors, of 

<ru = (2<r 2 + a f 2 )i (2) 

The relative motion between two features labeled 1 and 2 
is used to estimate velocity gradients and strain; we refer 
to it generically as stretching. It has the the form 

An s= u 2 - tii ( 3 ) 

and has an error with standard deviation 

aAu = (4ct 2 + 2 otf (4) 

If the points are close together, the geolocation errors for 
points 1 and 2 tend to cancel, and this error is reduced 
towards 

(2a f 2 )* (5) 

The values quoted for these errors have been: o g from 1 to 
5 km, and a from 0.9 to 2.5 km; see Table 2. A VHRR ice 
tracking has been performed almost exclusively with 
HRPT and LAC data with a resolution varying from 
1.1 km at nadir to 4 km at the swath edge. These values 
give a displacement error <r u (from Equation 2) between 
1.8 and 7.5 km, an upper bound for the stretching error 
oAu (from Equation 4) between 2.5 and 10.6 km, and a 

Table 2. Values for geolocation and feature 
identification errors , in km 


Reference cr g Reference <?( 


Clark and LaViolette 2 to 4 

(1981) 

Emery and others 1.5 Emery and others 0.9 
(1989) (1991b) 

Bordes and others 0.8-1. 7 Baldwin and Emery 1 

(1992) (1993) 

Heacock and others 2.7 Heacock and others 2.5 

(1993) ( 1993 ) 


lower bound for the stretching error aAu (from Equation 
5) between 1.3 and 3.5 km. These errors are quite 
significant compared to typical daily ice motions 
between several km and 10 or 20 km, and stretching of 
roughly 1 to 10km per 100km of spatial separation. 
Many of the low estimates of error magnitudes are for 
ideal cases: ice near land and near swath centers with a 
nadir view and the smallest image pixels. 

For AVHRR ice tracking in particular, it seems 
possible to make a clear statement about desired error 
size. If the geolocation and feature identification errors 
can each be limited to about 1 km, the motion data are 
quite useful scientifically. Measurements of this accuracy 
produce ice displacements with errors of 1.7 km and 
stretching errors (over a short baseline, from Equation 5) 
of 1 .4 km. 

7.5. Recommendations 

Because errors can render ice motion measurements 
scientifically useless, effort should be directed at docu- 
menting and testing errors more thoroughly, and in 
reducing them to their lower limits. It would be good 
practice to attach an error estimate to every ice motion 
vector and to treat this error as a function both of cross- 
track position within the swath, and of distance from 
ground points used to geolocate the image. 

1. The feature identification error <Tf increases with 
scan angle. Near the swath edges cross-track pixel 
size is 4 km. The practice has been to use data 
within a vaguely defined “center of the swath”, and 
ignore this variation. We recommend that feature 
identification error be treated as a function of cross- 
track position in LAC and HRPT data, and that 
this dependency be documented. 

2. Geolocation error a g is believed to increase with 
distance from ground points. Some test of this 
notion is needed. Tests could be run with con- 
tinental scenes, treating them as sea-ice images with 
land in only a portion of the scene. If errors are 
clearly shown to grow with distance from land, this 
behavior should be incorporated in error estimates 
for all ice motion vectors. 

3. The errors in ice motion data are generally errors in 
position, not errors in the time of the images. 
Dividing displacements by the time interval 
between images to compute velocity can obfuscate 
the discussion of errors. We recommend using 
position and displacement, not velocity, when 
discussing errors. 

4. Validation of the combination of geolocation and 
feature identification error should be accomplished 
by comparison with buoy drift data and with SAR 
satellite ice tracking. 


8. ICE SHEET STUDIES 

8.1. Introduction 

AVHRR data can be used for ice sheet studies for three 
main purposes: (1) surface temperature, discussed in 


Section 3; (2) surface morphology of the ice sheet; (3) 
calving of the ice fronts. Recently work has been done on 
deriving accumulation rate and snow grain-size from 
AVHRR data, but error estimations are needed before 
applying these methods over large areas of the ice sheet. 

8.2. Surface morphology 

Morphological studies include the identification of surface 
features (Casassa and others, 1991) and the characteriza- 
tion of different modes of ice flow (Bindschadler and 
Vornberger, 1990; Seko and others, 1993). This inform- 
ation can indicate present and past flow conditions of the 
ice and is especially important for improving the map 
coverage of large ice masses. 

Since surface features do not change rapidly with 
time, frequent coverage is not required, but accurate 
georeferencing is essential for determining the location. 
The fine spatial scale of these features, however, requires 
AVHRR of the highest spatial resolution (LAC near the 
swath center) and highest radiometric resolution (10-bit). 
With these resolutions, topographic roughness of approxi- 
mately 10~ 3 over horizontal scales of more than a few 
kilometers can be detected. Topographic features are 
revealed most effectively at orthogonal sun illumination. 
Thus, the need for imagery of the same area at various 
sun azimuths is recommended. The permanence of these 
surface features affords the opportunity for mosaicking 
different images to obtain the best surface information. 
Visible and near-infrared bands are best for this purpose. 

8.3. Calving Rate 

Monitoring ice fronts with AVHRR data on a periodic 
basis permits the detection of large iceberg calving events 
(Ferrigno and Gould, 1987; Keys and others, 1990). 
AVHRR data is also useful for monitoring the position of 
the ice front and comparing with historical data for 
deriving calving rates with accuracy of 0.2 km a -1 . With 
this kind of accuracy it is possible to monitor fast-moving 
ice fronts with AVHRR data over periods of 5-10 years to 
obtain reliable calving rate information. 

High spatial resolution is needed and accurate 
georeferencing is essential. Visible and near-infrared 
channels are preferred, although ice fronts also appear 
clearly in thermal infrared bands, which permits 
monitoring during winter at high latitudes. 

8.4. Potential for other ice sheet studies 

In the interior of East Antarctica AVHRR brightness 
patterns correlates strongly with surface topography 
(Seko and others, 1993). Since accumulation rate is 
closely related to topography, AVHRR data can be used 
to extend the sparse data on accumulation rate from East 
Antarctica. However, an assessment of the accuracy of 
this method is needed. The multispectral dependence of 
reflectivity on snow grain-size has permitted the quanti- 
fication of grain-size from Landsat-TM data (Bordelles 
and Fily, 1993). This method could be used with 
AVHRR visible and near-infrared bands for deriving 
snow grain-size spatially. This can provide an indepen- 
dent measure of surface nnA in 


Steffen and others: Snow and ice application of AVHRR 

extraction of other parameters (accumulation rate, 

physical temperature) from passive microwave data. 

Error estimation is needed for this method as well. 

8.5. Recommendations 

1. AVHRR of the highest spatial resolution (LAC 
near the swath center) and highest radiometric 
resolution (10-bit) are required to resolve surface 
features. 

2. Accurate georeferencing is essential for determining 
the location as ice sheets have only a few fixed 
reference points. 

3. Imagery of the same area at various sun azimuths 
and orthogonal sun illumination are needed to 
reveal topographic features most effectively. 

9. CONCLUSIONS 

The previous seven sections discussed the great potential 
AVHRR applications, but also showed major short- 
comings which need to be taken into account for polar 
applications. Today, a major emphasis in polar research is 
directed towards monitoring surface parameters that 
describe present climate conditions and help us under- 
stand how the climate system functions. AVHRR data 
have the potential to provide a rich data stream bearing 
on polar conditions and processes, but these data have 
been underutilized, due mainly to the inaccessibility of 
calibrated, geolocated LAC and HRPT data sets. Many 
applications, such as cloud masking, surface temperature, 
lead statistics, ice concentration, and ice motion, require 
multi-temporal data sets for process studies and for 
reducing the uncertainty in derived surface parameters. 
Initial steps to provide such data sets have been made by 
the ONR-supported AVHRR 1989 data set (Fetterer and 
Hawkins, 1993). 

Retrieval of ice surface temperature shows promise 
and would represent a significant contribution to the 
field, considering the current lack of spatially detailed 
measurements in the polar regions. This might be the 
most important application of AVHRR (and similar) 
data in the polar regions, since some of the other 
applications can be more effectively met by other satellite 
systems. But retrieved data have not been yet validated, 
although the estimated accuracy is believed to be between 
within 1 and 4K. Atmospheric models (radiative trans- 
fer) to account for absorption and emission of thermal 
radiation are in place, but additional information on ice 
crystal precipitation and aerosol distribution are needed 
to reduce the uncertainty in 1ST retrieval. The same 
parameters are needed for an improved albedo retrieval. 
Currently, narrow-band albedo can be retrieved for cloud 
free pixels within an accuracy of 1^1% for summer 
conditions (no ice crystal precipitation) if coincident 
atmospheric profile measurements are available. Ice 
concentration algorithms make use of surface reflectance 
and surface temperature information and consequently 
also need improved atmospheric input data to character- 
ize temporal and spatial variability of water vapor, 
diamond dust, and ice crystals in conjunction with 
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be placed on use of a combined sensor algorithm (i.e. 
passive microwave, thermal and visible, synthetic 
aperture radar) to circumvent present problems with 
cloud masking in the visible and thermal spectra. 

Lead detection has been highly successful for some 
individual case studies, but estimates of lead statistics and 
other parameters required as input to models have to be 
refined. Additional research is needed to understand the 
pattern of large-scale and small-scale leads and the 
variability in lead-width distribution (i.e., A VHRR in 
combination with Landsat or SAR data). Ice motion 
retrieval depends in the first order on the geolocation 
accuracy of the AVHRR data. This accuracy was 
reported to be between 0.8— 4km, depending on location 
and the presence of distinct geographical features. 
Feature identification is also crucial for ice tracking, 
and presently this accuracy is quoted to be between 0.9- 
2.5 km, also depending on the authors. Assuming a 
geolocation and feature identification error of 1 km, and 
ice displacement error of 1 .7 km, an ice stretching error of 
1.4 km could be expected. 

Having taken stock of AVHRR polar algorithms, we 
see three major objectives for advancing the scientific 
utility of AVHRR data. First, development of unsuper- 
vised cloud masking techniques and subsequent evalua- 
tions of their accuracy must be attacked with great 
energy. Cloud masking is essential for all retrievals of 
surface parameters, but today there is no operational 
polar cloud masking algorithm in place with sufficient 
accuracy (>90%) to be useful to the research and 
operations communities. Night-time cloud masking is far 
from being operational, and it is too early even to 
estimate the accuracy that may be achievable. In a real 
sense, the development of all-surface parameter algo- 
rithms for large scale studies is held up by the lack of 
adequate cloud masking. 

* We recommend that the different polar cloud 
masking groups work closely together to guarantee 
coordination and timely development of satisfac- 
tory algorithms. 

Second, clear statements of the accuracies of parameter 
retrievals must be provided by algorithm developers so 
that the scientific community can assess in the geophysical 
data derived from the sensor. Algorithms accepted by a 
knowledgeable portion of the community must be widely 
disseminated. The present situation is that the validity of 
some retrievals is open to question; the scientific 
community will not support devoting resources to the 
processing of data in which they have no basis for 
confidence. 

* We recommend that funding agencies devote more 
resources to validation work and to careful error 
estimation as part of research on these algorithms. 

Third, acquisition and processing of polar AVHRR data 
must become a coordinated efTort between the science 
community, receiving stations and data centers, and the 
funding agencies to guarantee a long-term commitment 

to the collection and distribution of high quality data. 
. • : associated 


archives should begin a dialogue among themselves and 
with their scientific users and funding agencies to establish 
uniform data standards and even common formats for 
geophysical data products. Algorithms to be used must be 
chosen by or with the consent of the scientific community. 
Such action will free geoscientists from the need to post- 
process satellite data, and allow them to concentrated on 
the study of geophysics. If such standardization among 
stations and archives is not possible, then a single site for 
archiving polar high-resolution AVHRR data should be 
identified and funded to collect data from various stations 
and to convert it to geophysical products of common 
formats and known accuracies. 

* We recommend a workshop be convened to provide 

a forum for representatives of all polar receiving 
stations and representatives of the scientific user 
community to investigate ways to implement a 
common data catalog, and to standardize AVHRR 
data products. 

* We recommend that a data center be identified to 
coordinate acquisition, processing, and distribution 
between all receiving stations. It would be desirable 
to have the current data acquisitions accompanied 
by suitable ancillary data sets that could be used to 
improve interpretation of the data and to validate 
derived products. 

* We recommend careful evaluation of information 
content of AVHRR LAC versus AVHRR GAC 
data. GAC data is more suitable for global analysis 
in the context of monitoring and climate change 
studies. 

Until such steps are taken, this valuable data source will 

continue in its present state of underutilization. 
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